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Abstract 

A dedicated test of the effects of Nitrogen contamination in liquid Argon has 
been performed at the INFN-Gran Sasso Laboratory (LNGS, Italy) within the 
WArP R&D program. 

A detector has been designed and assembled for this specific task and connected 
to a system for the injection of controlled amounts of gaseous Nitrogen into the 
liquid Argon. Purpose of the test is to detect the reduction of the Ar scintillation 
light emission as a function of the amount of the Nitrogen contaminant injected in 
the Argon volume. A wide concentration range, spanning from ~10~^ ppm up to 
~10^ ppm, has been explored. 

Measurements have been done with electrons in the energy range of minimum 
ionizing particles (7-conversion from radioactive sources). Source spectra at 
different Nitrogen contaminations are analyzed, showing sensitive reduction of the 
scintillation yield at increasing concentrations. 

The rate constant of the light quenching process induced by Nitrogen in liquid 
Ar has been found to be A;(N2)=0.11 /xs^^ppm"-^. 

Direct PMT signals acquisition at high time resolution by fast Waveform 
recording allowed to extract with high precision the main characteristics of the 
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scintillation light emission in pure and contaminated LAr. In particular, the 
decreasing behavior in lifetime and relative amplitude of the slow component is 
found to be appreciable from 0(1 ppm) of Nitrogen concentrations. 

Key words: Scintillation Detectors, Liquid Noble Gases, Quenching (fluorescence). 
Dark Matter Search 

PACS: 29.40.Mc, 61.25.Bi, 33.50.Hv, 95.35. +d 
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1 Introduction 

Dark Matter search in the form of WIMPs is of primary interest in the 
present astroparticle physics scenario. Direct detection of Dark Matter with 
noble gases liquified as target medium is one of the most promising line of 

* Corresponding authors 
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development in experimental technology. Argon is an ideal medium and the 
feasibility of Ar-based detectors has been firmly proved by the R&D study of 
the WArP Collaboration [IJ. The two-phase (liquid-gas) technology developed 
by WArP is based on the simultaneous detection of both signals produced by 
ionization events in liquid Argon (LAr): free electron charge and scintillation 
light. 

In the case of the ionization charge, the main limitation to its full collection 
comes from electron attachment process due to contaminations of electro- 
negative impurities like Oxygen present at residual level in (commercial) 
Argon. Detailed studies have been performed and led, since the pioneering 
work of the Icarus Collaboration [2], to the development of adequate 
purification systems (Oxygen reactants and molecular sieves, down to < 0.1 
ppb of Oxygen equivalent concentration), normally employed in experimental 
applications. 

On the other hand, the effects of impurities on the scintillation light yield 
are rather less precisely explored. Quenching (i.e. non-radiative) process in 
two-body collision of impurity-molecules with Ar2 excimer states, otherwise 
decaying with VUV light emission, may take place depending on the type 
of impurity and the concentration level. This led to include within the 
WArP R&D program investigations on Oxygen and Nitrogen contamination 
effects on LAr scintillation light. Two experimental tests have been thus 
performed. 

The case of Nitrogen (N2), present in commercial Argon at the level of 1 to 
10 ppm, depending on the selection grade, is the subject of the present study 
(while results from the test with Oxygen have been reported separately [3]). 
Quenching of the light yield in N2-contaminated LAr has been found, 
in particular of the slow component of the scintillation emission. This 
effect is relevant not only because it imposes some limitations on the full 
collection of the light available, but also because it deteriorates the detector 
capability of Pulse Shape discrimination of background events (electrons from 
7 interactions) from the signal (Ar-recoils, potentially induced by WIMP 
interactions) [T]. 



2 Features of the scintillation radiation in LAr 

It is experimentally established and theoretically well understood that 
recombination and de-excitation processes following the passage of ionizing 
particles in liquid Argon (LAr) finally originate scintillation radiation [1], [5]. 
In particular, interactions of ionizing particles in LAr cause the formation of 
both electron-hole (Ar+) pairs and Ar* excited atoms (the ratio Ar*/Ar"'" being 
0.21 [U]). Ar* excited atoms lead to the formation of the Arj low excited dimer 
through collision with Ar atoms ("self-trapping" process). On the other hand. 
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Ar"*" ions also lead to the formation of through a number of subsequent 
processes including recombination f7]. The excited dimer states formed in both 
cases in LAr are recognized to be the singlet ^Tj^ and the triplet excimer 
states in the M-band, characteristic of the solid Argon structure [Hj . The rise- 
time corresponding to excimer formation and relaxation is very fast, in the 
sub-nanosecond range for both components. 

The de-excitation processes to the dissociative ground state ^E^ lead to 
scintillation light emission in the Vacuum Ultra- Violet (VUV) region: 

Ar;^2 Ar + l-f (1) 

In liquid phase, various electron-to-ion recombination mechanisms along the 
ionization track are active, depending on the type of the ionizing particle and 
on its LET (Linear Energy Transfer, i.e. the specific energy loss along the path) 
[0]. These mechanisms affect the number of excited dimers Arg produced per 
unit of deposited energy, as well as the relative populations of the singlet and 
triplet states. For a minimum ionizing particle (mip) the photon yield has been 
measured to be 4.0 x 10"^ 7/MeV |6]. 

The 7 decay spectrum of both excimer states has been extensively investigated. 
It is well represented by a gaussian shape, peaking at A = 127 nm with 
FWHM ~ 6 nm (at boiling point T=87.3 K, p = L395 g/cm^) [9]. The 
time dependence of the photon emission is known less precisely. The various 
measured values of both the excimer lifetimes and of the relative amplitudes 
reported in literature are in fact quite different with each other. To a first 
approximation, in all measurements the overall scintillation light emission 
exhibits a double exponential decay form (see [9] for a complete compilation 
of available data), characterized by two very different components: a fast 
component, with a time constant ts in the 2 ns to 6 ns range, and a slow 
component, with a time constant tt in the 1100 ns to 1600 ns range. These 
are associated to the lifetimes of the singlet and of the triplet ^S^j states 
in LAr, respectively. 

In addition to these fast and slow components, an intermediate component 
(with decay time around 40 ns) has been also sometime reported in literature 
[TT] [T^ . the origin of which was never investigated in details nor definitively 
confirmed. 

While the time constants do not depend appreciably on the ionization density, 
the amplitude ratio As /At of the singlet and triplet states is found to be 
strongly affected by the ionization density pO]. In particular, all measurements 
show an enhancement of the formation {fast component) at higher 

deposited energy densit}p] As an example, the relative amplitude for the fast 
and for the slow component in case of a mip is reported to be As/ At = 0.3 [TT] 
(i.e. As = 23% and At = 77% respectively), while for heavily ionizing 

^ This effect is contrary to the effect in organic scintillators, in which the relative 
intensity of the slow component increases with increasing specific ionization density. 
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particles the intensity ratio increases (e.g. As/ At = 1.3 for a-particles and 
=3 for nuclear recoils, but higher values are given elsewhere [12] )• This wide 
separation is an important feature of the scintillation signals in LAr, leading 
to define robust Pulse Shape Discrimination criteria suitable for particle 
identification. 

The lack of accuracy/consistency in the presently available data can 
presumably be ascribed to (i) light quenching from residual concentration 
of impurities diluted in LAr, (ii) limited time resolution employed in past 
experiments, (iii) dependence on the LAr density (not systematically reported 
or accounted for). 



2. 1 Light quenching from diluted N2 impurities in LAr 



The de-excitation process rate can be described by a first-order rate law, 
characterized by two decay time constants ts and tt associated to the singlet 
and the triplet '^S„ excimer states of Ai^, respectively (j = S,T): 

= - I, [Ar;h [ArUt] = [Ar*],(0) e"*/-^' (2) 

Assuming that one scintillation photon per Ar2 dissociation is produced, the 
time dependence of the scintillation light emission in pure LAr can thus be 
represented by the following probability distribution function (pdf): 

m = — exp(--) + ^ exp(--) (3) 

Ts Ts Tt Tt 

with the sum of the relative amplitudes constrained to unity {As + At = 1) 
so that 

roo 

i{t)dt = {As +At) = I. (4) 







Light quenching due to (residual and unknown) impurity concentration in LAr 
is often indicated as the dominant source of uncertainty in the available data 
on LAr intrinsic properties. 

Residual concentration at the ppm (part per million) level of N2, O2, H2O and 
CO-I-CO2 contaminants (whose presence is usually indicated in commercially 
available Argon) can lead to a substantial reduction of the scintillation light 
intensity, mainly due to the decrease of the slow component amplitude. 
Argon purification systems (Oxygen reactants and molecular sieves) are known 
to be very effective in reducing the O2 and H2O, CO+CO2 contamination up 
to a negligible level (< 0.1 ppb, part per billion). At present, methods for 
removing the N2 residual content in Argon are instead less commonly used 
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in experimental application^^ and the effects of N2 contamination are rather 
poorly known. 

The main reaction to be considered is the quenching process in two-body 
collision of N2 impurities with excimer states [2] : 



Arl + N2 ^ 2Ar + N2 (5) 

This non-radiative collisional reaction is in competition with the the de- 
excitation process leading to VUV light emission. As a result, a sensitive 
quenching of the scintillation light yield is expected, in particular for lightly 
ionizing particles. In fact, the quenching process leads to the decrease of the 
excimer concentration [Ar2], while the contaminant concentration [N2] stays 
constant in time. To a first approximation, also for this case a first-order rate 
law can be assumed and characterized by the rate constant k: 

^^i^ = - km[Arl] [Ar*m = [Ar;](0) e-'^t^^l * (6) 

The value reported in literature is A; ~ 3.8x10^^^ cm"^ s~^ [H] (equivalent to 
0.11 ppm~^ fis~^). 

The time dependence of the scintillation light emission in N2-contaminated 
LAr can be represented by an expression formally identical to EqjS] 

nt) = ^ exp(-A) + ^ exp(--^) (7) 

where rj [j = S,T], the effective (decreased) lifetime of each component, is a 
function of the [N2] concentration: 

km) = 7 + ^ m (8) 

and A'j, the effective (quenched) amplitude, is also function of the [N2] 
concentration: 

Notice that under these definitions £'(t) reduces exactly to l{t) for [N2]=0, 
namely for pure Argon. The sum of the quenched amplitudes of Eq. ([T]) is no 
longer constrained to unity: 

;>oo 

/ e{t)dt = (A's +A't) < 1 (10) 
Jo 



^ O2 and CO+CO2 are known to be electro-negative molecules. Their high e-affinity 
induces ionization electron attachment, with dramatic effects on detectors based 
on charge collection (LAr-TPC). N2 molecules are not electro-negative [15], and 
therefore no appreciable effects on the free electron charge can be expected. 
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From Eq. (10) it can be defined the overall quenching factor (Qf)'- 



Qf = A's + A't ; < Qi. < 1 (11) 

namely, the ratio between the total intensity of scintillation light emitted for 
a given [N2] contamination with respect to the case of pure liquid Argon. Q f 
also represents by definition the normalization value for i'{t). 

This work is dedicated to a systematic study of the the effects of Nitrogen 
contamination in LAr. In particular, the aim is: 

• to determine the overall quenching factor affecting the LAr scintillation light 
yield at various N2 concentrations [Sect, [ij; 

• to determine the effective lifetimes rj as a function of the N2 contamination 
for each component, as well as the quenched amplitudes A'^ [Sect, [sj; 

• to extract the main characteristics of the scintillation light emission in pure 
LAr and the value of the rate constant k of the quenching process associated 
to the presence of N2 contaminant [Sectj5.1|. 



3 Detector and experimental method 



A detector has been designed and assembled at the INFN-LNGS external 
facility ( "Hall di Montaggio" ) for this specific task. It consists of a cylindrical 
LAr cell in PTFE (h=12 cm, 0=8.5 cm internal dimensions, wall thickness 
0.5 cm) containing about 0.7 It of LAr (about 1 kg of active mass), viewed 
by a single 2" photo-multiplier (Electron Tubes ETL D745UA) mounted on the 
top end of the cell. 

The photo-multiplier (PMT) is manufactured to work at cryogenics 
temperatures (LAr); however, the glass window is not transparent to the 
127 nm LAr scintillation light. Therefore, the PMT optical window is covered 
by a layer of TetraPhenyl-Butadiene (TPB), an efficient frequency down- 
converter with emission spectrum peaked at 430 nm and width of about 50 
nm (FWHM), well matching the PMT sensitivity range. To enhance the light 
collection a reflector layer (also covered with TPB) surrounds the internal 
walls (side and bottom) of the PTFE cell. 

The detector is housed in a stainless steel cylindrical chamber (h=85 cm, 
0=10 cm), closed at both ends by vacuum-tight ConFlat (CFIOO) flanges. 
Fig. [1] (schematic layout). The internal volume of the chamber is about 6.5 It 
and it contains, after fllling, a total amount of 3.0 It of LAr (including the LAr 
cell active volume at its bottom). The chamber is immersed in a LAr bath of a 
stainless steel open dewar, to liquify and keep at stable temperature the LAr 
internal volume. 

A transfer line for the Ar fllling and for the injection of controlled amounts 
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Fig. 1. Schematic layout of the experimental set-up for the A/2 contamination test. 

of gaseous Nitrogen (GN2) has been also assembled and connected to the 
LAr chamber via vacuum-tight pipes and feed-throughs. The line includes a 
two-stages system formed by a gaseous Nitrogen (GN2) buffer and a mixer of 
gaseous Nitrogen and Argon, connected via UHV valves, Fig. [Tj The internal 
volume of the buffer is precisely known and accurate pressure control allows 
the injection of controlled amounts of N2 gas in the LAr volume (from 1 ppm 
to 300 ppm of N2 per injection with < 6% precision). 

Gaseous Argon (GAr) for filling, supplied by a 200 atm pressurized bottle, is 
the best grade 6.0 (99.9999%) commercial Argon, with impurity concentration 
below 1 ppm. In Tab jl] the purity specification for the GAr in use is reported 



The GAr from the 6.0-bottle is flushed during filling through a 
Hydrosorb/Oxysorb cartridge positioned along the GAr line for partial 
removaf^ of O2 and H2O. The GN2 for the controlled contamination is 5.5 
GN2 type (99.9995%), also supplied by a pressurized bottle. 
In Fig. [2] a picture taken during the test is shown. 

After vacuum pumping of the whole system (detector chamber and transfer 



^ No recirculation system is implemented for further O2 purification. Initial content 
of residual N2 in the 6.0 Ar is not filtered out. Additional contaminations from 
material out-gassing inside the chamber cannot be excluded. 
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Table 1 

Argon 6.0 (from Rivoira Supplier) purity specifications: maximum impurity levels. 



Oxygen (O2) < 0.2 ppm 

Nitrogen (N2) < 0.5 ppm 

Water (H2O) 0.5 ppm 

Hydrogen (H2) 0.1 ppm 
Total Hydrocarbons (THC) 0.05 ppm 

Carbon Dioxide (CO2) 0.05 ppm 

Carbon Monoxide (CO) 0.05 ppm 




9 



line) up to a residual vacuum of few 10~^ mbar, the external bath is filled with 
commercial LAr, the GAr line is opened and the purified Ar gas feeding the 
detector chamber is liquified inside it up to reaching a given level, precisely 
monitored by a level meter, corresponding to 3.0 ±1.5% It of LAr. 

At LAr working conditions of T=86.5 K, P=930 mbar (p^r=l-40 g/cm^), 
contamination with GN2 is performed by a controlled procedure: (1) GN2 
corresponding to a given fraction of the LAr molar quantitjj^ in the chamber 
is first transferred into the GN2 buffer volume (at controlled temperature and 
pressure), (2) the valve to the GN2+GAr mixer is opened, (3) GN2 is injected 
into the chamber by repeated fiushes of purified GAr through the mixer. 
Ar and N2 form almost ideal vapor/liquid mixture^^ In the liquid phase they 
are completely miscible and almost ideal liquid solution. The vapor pressure in 
the top part of the chamber (about 3.5 It above the liquid level) is different for 
Argon and Nitrogen and can be determined by the Antoine equation assuming 
the gas at the same temperature of the liquid Ar. Nitrogen results to be more 
volatile than Argon by a factor of 3; therefore the Nitrogen concentration 
in gas phase is higher than in the liquid [T7|. However, since the gas/liquid 
volume ratio for Ar is very large (e.g. 835 at STP), the fraction of N2 injected 
in the chamber which is diluted in the LAr volume is very close to one. 

3.1 Data taking 

The PMT read-out system is structured as an oscilloscope channel: in fact, the 
PMT+read-out chain can be represented by the simplified scheme of Fig. |3} 
where at the first stage the PMT signal is duplicated by a linear Fan- out, 
one output being used as input for a Constant Fraction Discriminator (CFD) 
providing a trigger signal when a defined threshold is passed, and the second 
output is directly recorded by a fast Waveform Recorder (Acqiris, DP235 Dual- 
Channel PCI Digitizer Card). At each trigger from the CFD the signal waveform 
is recorded with sampling time of 1 ns over a full record length of 10 /is. A 
LabView application has been developed for the data acquisition and storage. 
In Fig. |4] a typical waveform recorded during the N2 test is shown. 
After the filling procedure with purified 6.0 Argon was completed, the N2 
experimental test started on February 2007 and lasted about one month. The 
test was performed by adding progressively controlled amounts of N2 and 
exposing the LAr cell to 7-sources after each contamination. Data samples 
(10^ waveforms per run, trigger rate ~ 0.5 kHz) have been recorded in two 

^ The concentration units in use here are ppm, parts per million (atomic, often 
indicated as ppma). This indicates the ratio between the number of interesting 
elements (in our case N2 molecules) to ordinary elements (Ar atoms). 
^ The dilution process is rather slow. However, the speed is substantially increased 
by the repeated flushing of GN2+GAr mixture through the liquid. 
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LAr 
Cell 



Reflector Layer with TPB 



Fig. 3. PMT Electronics read-out chain. 




Time (ns^ 



Fig. 4. Typical (single) waveform recorded during the N2 test. Event with large energy 
deposition from cosmic muon (mip) crossing the LAr cell. 



independent runs per N2 concentration, with ^^''Cs and ^°Co 7-ray sources. 
After the first run at ppm (i.e. no additional N2 injection to the initial 
purified 6.0 Argon), the contamination levels were set at [N2] = 1 ppm, 2 
ppm, 7 ppm, 12 ppm, 20 ppm, 40 ppm, 60 ppm, 100 ppm, 500 ppm, 1000 ppm 
and 3000 ppm. 

Single photo-electron measurements were also acquired before and after each 
7-source run, to provide calibration data (see App. |A]) useful at various stages 
of the analysis performed. 
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4 Signal Amplitude Analysis 



Gamma rays from the ^°Co source (1.27 and 1.33 MeV) induce mainly single 
Compton interactions in the LAr cell active volume (Compton edge Ec — 1 
MeV), with final state electron in the mip range. Scintillation light is emitted 
and, after down-conversion and refiections at the active volume boundaries, is 
collected at the PMT photo-cathode; finally the signal waveform is recorded. 
By single waveform integration, the absolute (individual) signal amplitude is 




TOO BOO 9DD 

energy <phe) 

Fig. 5. ^^Co source Compton spectra from each [N2] run. 

obtained (in units of photo-electrons, eventually proportional to the electron 
energy deposited in the LAr cell). 

Pulse amplitude spectra have been thus obtained for each run at different 
[N2] value, as reported in Fig. [sjwith the ^"Co source. Manifestly, the spectra 
are progressively down-scaled at increasing N2 concentration as due to the 
quenching process of reaction 

A dedicated fitting procedure has been developed to determine the scale factor 



Qf (the overall Quenching Factor defined in Eq, 11) acting on the absolute 
signal amplitud^^ as a function of the N2 contamination. As an example, 
in Fig. [6] the result from the fitting procedure for the [N2]=20 ppm run is 
reported {Qf{[20 ppm]) =0.405: i.e. only ~40% of emit a photon, while 
the remaining ~60% undergo non-radiative N2 collisions). 



The Qf values and associated errors obtained from the analysis of the whole 

® Individual entry values of the uncontaminated (0 ppm) spectrum are down-scaled 
by Qfi free parameter, to best fit the experimental [N2] contaminated spectra. 
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Fig. 6. ^"Co source for the [N2]=20 ppm run (dashed histogram). Superimposed the 
result of the fitting procedure (red histogram). 

set of ^°Co runs at the various N2 contaminations are given in Fig|7l After 
a fast drop up to ~ 100 ppm {Qf{[100 ppm])=0.25) a slower decrease is 
observed at higher N2 concentrations (notice the log-log scale in use). The 
hght contribution from the slow component is expected to be vanished above 
[N2] > 100 ppm, and the fast component starts to be progressively affected 
by the quenching process because of the very high Nitrogen concentration. 

The errors on the plot are from the fit (statistical) and are appreciable at the 
highest concentrations. 



5 Signal Shape Analysis 

As reported in SecjSj data were collected by full signal waveform recording 
(1 ns sampling time over a time window of 10 /is) for each trigger event. 
This allows a detailed reconstruction of the signal shape, in particular of the 
individual components (relative amplitude and decay time) of the scintillation 
light following an ionization event in LAr. 

To this purpose the recorded waveforms (wfm) have been appropriately 
treatecQ and averaged {V{t) for each [N2] run, 10^ events) and the averaged 
waveform amplitude (area) normalized to the corresponding Qp value. 

^ Individual waveform (wfm) are off-line processed: (1) 0.5 /is of baseline presamples 
are used for baseline determination and subtraction, (2) cuts on wfm amplitude 
(integral) are applied to remove PMT saturations and low energy events, (3) wfm 
are time equalized at the same peak position and finally (4) selected wfm's are 
summed together. 
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Fig. 7. Quenching factor Qp as a function of the [N2] contamination (^^Co runs). The 
value corresponds to the fraction of Ar\ surviving A/2 quenching and producing VUV 
photons. 

The averaged signal can be symbolically expressed as: 

v{t) = s{t) ® 7^(^) (12) 

where is the Response Function of the read-out system 

(PMT+cable+wfm digitizer), i.e. the electronic impulse response defined by 
the output signal obtained from an input single photo-electron pulse (usually 
indicated as SER, Single Electron Response). 

The true light signal shape S{t) can be obtained from the recorded waveform 
V(t) by a standard deconvolution procedure [T3], using the SER function 7l(t) 
experimentally obtained (as explained in App. |A]). In Figjs] three averaged 
signals (at [N2]= ppm, 12 ppm and 40 ppm) after SER deconvolution are 
shown. 

The signal function i'{t) (see Eqj?]) convoluted with a gaussian function has 
been used to fit the signal shape S{t) from each [N2] run. The gaussian spread 
(free parameter 0"^ ~ 4 ns from fit) takes into account the time resolution of 
the light collection at the detector sensitive area, due to transit time of the 
PMT, wave length shifting time, time spread due light propagation including 
multiple refiections at the volume boundaries, etc. 

Attempts to perform a fit with a signal function characterized by two time 
components (singlet and triplet) failed, as evident from Figjo} This led to use a 
three components model (S, T and I, Intermediate) , nicely fitting the present 
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Fig. 8. Signal shape at ppm, 12 ppm and 40 ppm o/" A/2 contamination, superimposed 
fits with three components. 

data (color lines in Figjs]). The presence of an intermediate component is thus 
confirmed, in agreement with some earlier indications [TT]. Dedicated studies 
on the origin of it are under way. Based on preliminary considerations, reported 
in App. [B} possible instrumental effects appear as the main contribution to 
this component. 



In Fig, 10 the quenched amplitudes of the individual components A'^ from the 
fit are shown as function of N2 contamination. The total amplitude (sum of 



the A'j amplitudes at given [N2], black stars in Fig 10) is constrained to the 
corresponding Qp value (Figj?]), as due to the average waveform normalization. 
The decay times rj of the individual components from the fitting procedure 



of the signal shapes at different N2 concentrations are reported in Fig, 11 



5. 1 Characteristics of LAr scintillation light and of N2 quenching process 



The quenched amplitudes and the decreased lifetimes of the scintillation light 
in N2-contaminated LAr are correlated variables depending on the intrinsic 
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Fig. 9. Signal shape at ppm contamination (black histogram) with comparison of 
two-components (red line) vs. three-component (blue line) model. 

LAr scintillation characteristics and on the rate constant of the N2 quenching 



process, as reported in Sec 2.1 



An overall fit of the data obtained from the signal shape analysis reported in 
Fig, 10 and Fig, 11 has been performed using Eq|9]and EqjSjas fitting functions, 



respectively. The values of eight free parameters have been thus obtained by 
X^-minimization: the LAr scintillation characteristics (intrinsic lifetimes tj and 
relative amplitudes Aj, j = S,T,I for mip electronic recoils), the rate constant 
k of the quenching process due to N2 contaminant, and the initial Nitrogen 
concentration [N2]m- This last corresponds to the N2 content in the 6.0 grade 
Argon used for the initial filling. Its value is obtained by using an additive free 
parameter in the fitting functions: [N2]— >[N2] + [N2]m. 

A contamination of Oxygen at the level of [O2]=0.06 ppm has been taken into 
account into the fit, as determined by the analysis reported in [3] from the 
parallel test of Oxygen contamination effects. 

The values of the parameters obtained from the overall fit are reported in 
Tab.2. 



Table 2 

Results from the overall fit: characteristics of scintillation radiation in LAr from mip 
electronic recoil and of quenching process from Nitrogen contamination. Quoted errors 
are statistical from the fit. 



Lifetimes 


Amplitudes 


N2 rate constant 


Initial N2 concentration 


Ts = 4.9±0.2 ns 


As = 18.8% 






Ti = 34±3 ns 


Ai = 7.4% 


k = 0.11 ±0.01 ppm"^ iis-^ 


[N2]m = 0.40±0.20 ppm 


TT = 1260±10 ns 


At = 73.8% 
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The rate const anl| ^ |of the N2 quenching process {k value) is in agreement with 
an earher measurement found in hterature 

The initial N2 concentration is within the expected range from the 
specifications of the 6.0 grade Ar used for filling. 

The relative amplitude of the intermediate component is small {Aj = 7.4%) 
compared to the amplitude of the other two components. It is also worth 
noticing that the amplitude ratio between As + Aj and A^ is 0.35 in agreement 
with available reference data for light ionizing particles [TT] . 

In all components the fit is found to deviate from the data at the highest 
N2 concentrations, as possibly due to incomplete absorption of impurities 
into the liquid. This trend might be explained by a saturation effect of the 
solute (Nitrogen) in the LAr solvent. Similar indications have been reported 
by other groups [H] and also in our test with O2 contamination [3]. The 
measurements were performed on relatively short time scale (hours to day per 
contamination). This prevents to draw any definite conclusion on long term 
effects. 

A different parameterization depending on the N2 concentration under 
saturation hypothesis has been thus used for the fitting function of the lifetimes 
(each component): 

1 1 

- = - + k(3{l - e-^) (13) 

and of the amplitudes: 

4 = ^ (14) 

l + Tjk P {I - e~~) 

where the f3 parameter represents the concentration scale where saturation 
becomes effective [jS = 530 ppm from fit). At first order approximation both 
formulas reduce to the original fitting functions (Eqjsjand Eqjo]). 



The curves in Fig, 10 and Fig, 11 (full line in corresponding color) show the 
result of the fit with the modified parameterization. This model improves the 
scaling with N2 concentration and decreases the overall without changing 
at all the results reported in Tab. 2. 



^ The rate constant of the two-body (N2-Ar2) process is assumed to be the same 
for both the singlet and the triplet Ar excimer states. A fit has been tried 
with two rate constants as free parameters. This indicates that in case of N2-^Sm 
couphng, a slightly higher constant rate is preferred. 
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Fig. 10. Effective amplitudes (mip electronic recoil) of the individual components 
= S,T,I in corresponding color an symbol) as function of A/2 concentration. 
The total amplitude (sum of the A'j amplitudes at given [N2]) is constrained to the 
corresponding Qp value (black stars) obtained in Secj4[ Lines in corresponding color 
are from the fit with saturation model (Eq^^. 
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Fig. 11. Effective decay times of the individual components = S,T,I in 

corresponding color an symbol) as function of N2 concentration. For [N2]>500 ppm, 
^ (not shown). Lines in corresponding color are from the fit with saturation model 
(Eq\l^. 



18 



The errors associated to the results in Tab. 2 are from the fit (statistical). 
These take into account the statistical errors from the fit of the waveforms 
at different N2 contaminations and also the error associated to the injected 
amounts of Nitrogen in the contamination procedure. 

The goodness of the overall fit (saturation model) is satisfactory (C.L.~80%). 
Systematic uncertainty may come from different sources: instability of the 
signal read-out system, bias due to residual concentration of unknown types of 
impurities in LAr leading to light absorption or further Ar-excimer quenching, 
possible uncertainty introduced by the procedures of data analysis. Only this 
last issue has been investigated in some details. 

The off-line analysis of the light signal shape (obtained from the PMT 
waveform recording at the different N2 contaminations) is based on the 
deconvolution procedure of the SER function (see AppjA]). Systematic 
uncertainty associated to this procedure is relevant mainly to the long- 
lived component determination and has been estimated to be around 4%: 
[tt = 1260 ± 10 (stat) ± 55 (syst) ns]. Comparisons with alternative fit 
procedures confirm this evaluation, as reported in details in AppjCj 



6 Conclusions 

New generation detectors based on LAr as active medium for Dark Matter 
direct search and for neutrino physics exploit both scintillation light and free 
electron charge from ionization events. 

Residual content of electro-negative impurities in LAr (like Oxygen molecules) 
significantly reduces the amount of charge and light by attachment and 
quenching processes respectively. This led to the development of dedicated 
O2 purification systems, routinely employed in the experimental set-up's. 
On the other hand, non electro-negative contaminants like Nitrogen can also 
be found at appreciable concentration level even in best grade commercial 
Argon. Contamination at ppm level of N2 leads to a substantial reduction 
of the scintillation light intensity, due to the quenching process in two-body 
collision of N2 impurities with Ar2 excimer states otherwise decaying with 
VUV light emission. (No appreciable effects from N2 contamination on the 
free electron charge can be expected, due to the low electron affinity of the N2 
molecule). 

A dedicated test has been performed by means of a controlled N2 
contamination procedure. Measurements have been done by mip particle 
excitation (Compton electrons from 7-sources in the MeV range). 
The effect on the scintillation light collection has been measured over a wide 
range of N2 concentration, spanning from ~10~^ ppm up to ~10^ PPm, 
though a saturation effect of the solute (Nitrogen) in the LAr solvent has 
been presumably found at the highest contaminations. 
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The rate constant of the hght quenching process induced by Nitrogen in LAr 
has been found to be A';(N2)=0.11 /xs^^ppm^^ (in agreement with some early 
measurements reported in hterature). This imphes that, for example, at ~ 1 
ppm level a ~ 20% reduction of the scintillation light is experienced due to 
the N2 quenching process. The quenching factor affecting the light yield has 
been determined over the whole explored range of N2 concentrations. 

The experimental test also allowed to extract the main characteristics of the 
scintillation hght emission in pure LAr, thanks to the direct PMT signal 
acquisition by the fast Waveform Recorder in use with the implemented DAQ 
system. 

The light signal shape is well represented by the superposition of three 
components with exponential decay. The fast and the slow components, 
with decay time constants T5=4.9 ns and tt=1260 ns, are recognized to be 
associated to the decay of the singlet and the triplet excimer states 
of Ar2. An intermediate component is also found to be present (as sometime 
reported in literature), whose origin could presumably be ascribed to PMT 
instrumental effects. 

The main effect of residual N2 in LAr is of reducing the slow component 
lifetime at increasing concentration, and consequently of varying the ratio of 
the relative amplitudes to higher values. This makes the slow signals from 7/e 
background less distinguishable in Pulse Shape vs the fast Ar-recoils signals 
possibly induced by WIMP interactions, in Dark Matter search experiments. 
Therefore, the implementation of dedicated methods for removal of the 
residual N2 content results to be recommended with LAr based detectors, 
or at least the use of best grade commercial Argon with reduced Nitrogen 
contamination (below ~0.5 ppm) appears as definitively necessary. 
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A Single Electron Response 



Single photo-electron pulses (SER data) from thermionic dark counts have 
been routinely acquired during the test, before and after each source run at 
the different N2 concentrations. 
Pulse amplitude spectra from SER waveform integration were used for 

» 900 r 1 
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total charge (A. U.) 



Fig. A.l. SER spectrum of the PMT in use and fit superimposed. The interval around 
the mean value for single photo-electron waveform selection is also indicated. 

calibration purposes. Typical SER spectrum of the PMT in use is reported 



in Fig, A.l the Gaussian distribution at higher ADC values corresponds to 
the genuine single photo-electron mean amplitude and spread. The mean 
value (gaussian fit) gives the calibration constant, from ADC units into 
photoelectrons units, useful for Compton spectra analysis (Secjij) and checking 
of the detector stability in time. 



The single photo-electron pulses collected by full signal waveform recording 
(1 Gsample/s over 10 fis time window) were used for the determination of the 
read-out electronics Response Function: the impulse response TZ{t) is obtained 
by processing and averaging the SER waveforms with amplitude in a defined 
interval around the mean value of the corresponding single photo-electron 
pulse amplitude spectrum. 

A typical SER waveform averaged on a large sample of recorded signals is 
shown in Fig A. 2 This is characterized by two main features: (1) the genuine 
single photo-electron shape from the PMT with FWHM ~ 10 ns and a decay 
tail extending up to about 150 ns, and underneath (2) a very long, slowly 
decreasing component up to ~ 3/is (visible in the up-right, log scale panel of 
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Fig. A. 2. Example of SER averaged waveform TZ{t). 



Fig, A. 2) followed by a flat baseline. 

The flat baseline is due to pile-up (in the average wfm) from the PMT dark 
counts randomly distributed in the recorded time window, while the origin 
of the long-decay component is not fully understood. It may be ascribed to 
additional pile-up of after-pulses from ionization of residual gases in the PMT 
(normally reported in PMT test reports [IE])- Typically, these type of after- 
pulses are delayed in a (9(0.5 - 5) /xs interval after the main signal from which 
they originate, in the ratio of < 5% to the number of true pulses. 

The SER function lZ{t) (including the pile- up components) determined 
for each source run is used in the deconvolution procedure to obtain the 
corresponding true hght signal shape S{t) (Secj5|. 

The deconvolution of the response function 7^(t) from the experimental 
waveform V{t) is performed numerically in the frequency domain. The discrete 
Fourier transforms Vi and fj of Vi and IZi respectively are first evaluated by the 
Fast Fourier Transform (FFT) algorithm [T3]. According to the convolution 
theorem, the true light signal Si is hence obtained by taking (with the same 
FFT algorithm) the inverse Fourier transform of Si (= Vi/fi). 



B The intermediate component: preliminary considerations 



As described in SecjS} evidence for an intermediate component in the 
scintillation signal shape has been found at all the N2 contaminations, 
including the "0 ppm run" with lifetime tj ~34 ns and relative amplitude 
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Aj ~7.4%. The origin of this component is controversiah it could be 
ascribed to (1) instrumental effects (in the PMT), (2) presence of residual 
N2 concentration in LAr, (3) physical (unknown) properties of Argon. 
Explanations in terms of: 

- Intrinsic physical properties of Argon yet unknown or only marginally 
investigated cannot be excluded. For example, from the delayed decay of 
vibrationally hot Arg molecules surviving vibrational relaxation [8]. This 
decay, seen in gas at low pressure, was never investigated in liquid phase. 

- Light emission from Nitrogen trace in LAr is possible though rather unlike 
in the present experimental conditions: (a) scintillation emission from excited 
Nitrogen molecular states N2(C^nu) formed by energy transfer interaction 
with Ar* atoms, as seen in gaseous Ar-N2 mixtures, should be negligible since 
in liquid phase the precursor Ar* atoms are quickly self-trappecQ leading to 
Ar2 formation (moreover no increasing Aj amplitude is found at increasing N2 
concentrations in Ar-N2 studies [this paper] and [H]), however (b) a series of 
mechanisms leading to the formation of excited Nitrogen atomic states from 
Ar2 energy transfer has been recently identified and detected in Ar-N2(few 
%) gas mixtures The N* radiative decay (149 nm and 174 nm) may 
thus populate the intermediate component amplitude. The size of this effect, 
presumably small in liquid mixtures with very low N2 concentrations, as in 
the present study, is of interest for further studies. 

- Instrumental effects to explain the onset of the Intermediate component is 
our preferred option. Such effects could be ascribed to the use of PMT and 
enhanced by particular experimental conditions. For example, the so called 
anode glow in the PMT [18] can be considered: it has long been known that 
electrons in PMTs induce light emission in the last cascades of dynode systems. 
Photons may take the way back to the photo-cathode and so give rise to after- 
pulses, typically 20 to 50 ns after the true pulse from which they originate (and 
up to 120 ns in large PMTs). This effect has been well characterized in various 
studies, indicating that the total amount is of the order of 1%. 
Contributions from late-pulses due to photo-electron backscattering effect on 
the first dynodes can also be taken into account. Late-pulses are part of the 
main pulses but delayed by 5 to 80 ns, in the ratio of < 5% to the number of 
true pulses JH]- 

Both contributions may increase in high photo-electron density conditions, as 
in our case induced by the LAr fast decay component in the first few ns of 
scintillation emission. 

All these are preliminary considerations. Only dedicated studies (outside the 
limits of this report) may possibly indicate the true origin of this component. 



^ The Ar* self-trapping reaction rate is ~10^^ [H], while the collision frequency 
of the energy transfer process is ~10^ s~^ at [N2]=0.5 ppm (e.g. like in our "0 ppm 
run"). 
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C Systematic error in the slow component lifetime determination 



After completion of this report and while proof reading, an exhaustive paper on 
Scintillation time dependence and pulse shape discrimination from Lippincott 
et al. |2n] appeared in publication. In this paper a value for the slow component 
lifetime of rT=(1463 ± 5stat ± 50syst) ns is quoted. This value is only marginally 
compatible with the one we found here of rT=(1260 ± lOstat ± 55syst) ns. 
We ascribe this discrepancy to the different strategies adopted in the analysis 
of the LAr scintillation waveforms and not to physical or instrumental effects. 
In fact, by applying the procedure described in [201 to fit the tail of the average 
waveform V{t) obtained from the "0 ppm run" data sample we find a tt value 
for pure LAr in full agreement with the one reported in \2U\ . 
In more details, the value of the slow component lifetime we obtain with the 
fit model from [20] (characterized by an additional fiat baseline component) 

V{t) = Aexp{-t/TT) +B (C.l) 



is Tt =(1306 ± 5) ns. Taking into account that the Argon we used to fill 
the cell ("0 ppm run") had indeed some residual contamination of Nitrogen 
and Oxygen (400 ppb of N2 and 60 ppb of O2) and having determined the 
values of the quenching rate constants k for N2 [this report] and fc' for O2 [3] 
respectively, we can extrapolate the tt value for pure LAr, that results to be 
rr=(1453 ± 10) ns. 
In Fig, C.l we show the average waveform V{t) with the fit (C.l) 



Average Waveform with Cl-model fit 

Signal Siiape (from SER deconvolution) 
witli exponential fit 




time (iii) 

Fig. C.l. Fit comparison: average waveform V{t), as built by summing up tine waveforms 
recorded in the "0 ppm run", with C.l-model fit (green line) and signal shape S{t) 
(standard deconvolution procedure) with exponential fit (red line). 
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superimposed. For comparison the light signal shape S(t) as obtained 
from V{t) with our standard procedure based on the SER function 7l{t) 
deconvolution is also shown. The difference in the slope of the long-component 
is quite evident. In our data this is due to the presence of an "extra- 
long tail" exponentially decaying with lifetime of ~3 fis (rather than the 



simple flat baseline component), also present in the SER waveform (Fig, A. 2). 
This possibly may fake a longer triplet state decay time in the average 
non-deconvoluted waveform. Thanks to the deconvolution procedure this 
component is washed out. 

It is worth noting that when alternatively a fit of V{t) with an extra-long decay 



component is performed [B ^ B exp{—t/TB) in Eq C.l, with ~ 3 /is) the 
triplet component lifetime is determined within 2% spread around the value 
obtained from our standard procedure (SER deconvolution). 

To further explore possible sources of systematic uncertainty due to the off- 
line procedures in our data analysis, we also tested the classic method of 
Coincidence Single Photo- electron Counting [2T] for another determination of 
the long component lifetime. 

The Coincidence Single Photo-electron Counting technique has been largely 
exploited in many reports on scintillation lifetime measurements (for LAr see 
[II]5[I2],[I1])- Data recorded with our set-up (one single PMT with direct 
signal acquisition by fast Waveform Recorder) allow to perform an artificial 
single photo-electron counting experiment, based on off-line data treatment. 
Starting at 0.5 /xs after the onset of a triggered signal, a single photo-electron 
identifying algorithm was run through the recorded waveform V{t) and for 
each photo-electron pulse (defined by appropriate cuts) found at time a 
single bin with the value of 1 was substituted. 



V(t) Average waveform 
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Fig. C.2. Result of the artihcial single photo-electron counting method for tt 
determination (red histogram) with the "0 ppm run" data sample. The average 
waveform V{t) (black histogram) is also shown for comparison. 
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As an example with the "0 ppm run" data sample, the resulting histogram, 
which is per se cleaned from noise contributions, provides a tt determination 
that is in very good agreement (within ~ 2%) with the result of our standard 



procedure based on the SER function deconvolution, as shown in Fig, C.2 
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